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INTRODUCTION 


Thi» report summarizes research activities conducted Cron February 
1, 1985 to July 31, 1985 and pr el la inary conclusions regarding research 
objectives. The overall objective is to determine the feasibility of 
using Landsat data t o estlaate effective hydraulic properties of soils. 
The general approach is to apply Eagleson's "climatic-climax" hypothesis 
(Eagleson, 1982) to natural water-limited vegetation systems using canopy 
cover estimated from Landsat data. 

Natural water-1 lmlte<i systems typically consist of inhomogeneous 
vegetation canopies interspersed with bare soils. The ground resolution 
associated with one pixel from Landsat MSS (or TM) data is generally 
greater than the scale of the plant canopy or canopy clusters. Thus a 
method for resolving percent canopy cover at a subpixel level must first 
be established before the Eagleson hypothesis is tested. 

Specific project objectives are thus, 

1. Evaluate the existing capability of defining vegetation canopy 
density from Landsat data with a resolution and accuracy sufficient for 
use in tenting the climatic-climax hypothesis of cllmate-soil-vegetatlon 
equilibrium. 

2. Design. the research needed to create the necessary capability if 
it is found (under Objective 1) to be deficient. 

3. Design the series of experiments necessary to test the 
climatic-climax hypothesis. 

One major conclusion of the work to date is that current methods 
for resolving mixed pixels, which have been developed primarily for 
agricultural lands, have not been adequately tested on natural landscapes. 


Method* which require the ldentlf lee t Ion of training samples can not be 
easily applied to naturally vegetated ays tee* due to the lnhomogeneity of 
the cover. 

Two foreulatlon* are proposed which extend existing Methods of 
analyzing mixed pixels to naturally vegetated landscapes. The first method 
Involves use of the normalized vegetation Index. The second approach Is a 
physical model based on radiative transfer principles. Both methods will 
be analyzed for their feasibility on selected sites during the second half 
of this research project. 

RESEARCH ACTIVITIES 

The reporting period covers the Initiation of work on this topic. 
Accordingly, research activities have Included a visit to the Goddard Space 
Flight Center, Greenbelt, Maryland, comprehensive literature review, the 
preliminary formulation of algorithms to determine percent canopy cover at 
a subpixel scale of resolution, and the search for several representative 
watersheds In which to conduct the research. Details of each of these 
activities are described below. 

Visit to Goddard Space Flight Center, Beltsvllle, Maryland 
Research commenced with a visit by the principal Investigator, 
Professor Peter S. Eagleson, and his assistant, Michael F. Jaslnskl, to the 
Earth Resources and Hydrological Sciences Branches, Goddard Space Flight 
Center (GSFC) on January 31, 1985. The purpose of the visit was to 
familiarize the MIT researchers with ongoing activities at GSFC and to 
discuss research goals with key GSFC personnel. Meetings were held with 
Dr. Vincent V. Salomonson, Dr. Robert Gurney, and Dr. Thomas J, Schraugge, 
among others. Image analysis of MSS data was demonstrated by James 


Or a* by. During the afternoon, Profasaor Eaglason presented a seminar on 
water-limited equilibrium of savanna vegetation systems, work supported 
under a prior NASA research grant. 

Literature Review 

Ongoing literature review has concentrated on assessing the current 
capability of estimating percent canopy cover using Landsat data. Curraat 
techniques for classifying agricultural lands including those over mixed- 
crop types have been reviewed. The reflectance properties of vegetation 
canopies and soils in the visible and near Infrared range have been 
Investigated. This has required a review of the theoretical principles of 
radiative transfer through leaves, the atmosphere, and other media. 
Vegetation Canopy and Soil Reflectance . 

The reflectance of the visible and near infra-red range Is well 
documented. Canopy reflectances are typically much less than those 
measured for Individual leaves (Dickinson, 1983) due to the scattering 
properties of plant architecture and soil background. Colwell (1975) 
discussed the Important parameters which influence canopy reflectance, 
Including leaf hemispherical transmittance, leaf area and orientation, 
plant structure components, soil reflectance, solar zenith angle, look 
angle, and azimuth angle. Tucker et al. (1977) found that the effect of 
underlying soil decreases as vegetation biomass Increases. The results of 
numerous studies of plant and canopy reflectance are summarized in 
Rosenberg et al. (1983), Smith (1983) and Myers (1983). 

Radiative transfer formulas for homogeneous canopies have been 
developed in terms of the physical and optical properties of the plant and 
background soil reflectance (Suits, 1972; for summary of other formulas see 


Smith, 1983). Typically, homogeneous canopies have been modeled as a 
diffusing medium with absorbing and scattering properties. Dickinson 
(1983) applied the two-stream approximation for radiation transfer in the 
atmosphere (Meador at al., 1980) to plant canopies employing the Leaf Area 
Index (LAI) as a measure of the optical depth. 

The reflectivity of plant canopies is highly wavelength dependent. 

It Is typically low in the visible spectrum (<30%) and high in the near 
Infrared range (>50%). Kondratyev (1969) and Iqbal (1983) provide 
summaries of reflectances of various natural vegetation covers. 

The spectral characteristics of bare soils have also been extensively 
Investigated, as summarised by Meyers (1983). Soil reflectance Is a func- 
tion of color, mineral content, organic matter, particle else, texture, 
moisture content, and wavelength. Soil reflectance Is generally less than 
30% and Increases linearly with wavelength over the range 0.4 - 1.0 pm. 

Wet soils exhibit lower reflectance than dry soils but generally maintain r. 
similar linear behavior in reflectance over that range (Tucker et al., 
(1977). 

Reflectance from Natural Landscapes 

The reflectance from natural landscapes results from the complex 
Interaction of light with the vegetation canopy and soil types. The canopy 
alone often consists of a mixture of various plant species of different 
sizes and architectures. Soil types can also vary over a short distance. 
The radiance observed by a Landsat satellite is a composite of those canopy 
and soil components which are subsequently scattered and attenuated by the 
atmosphere. 

The Identification of landscapes which consist of several cover 
types is a common problem encountered by Image analysts. Both empirical 
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class if lection schemes as well as aora physically-based aodals hava bsen 
developed, although aaplrlcal asthods ars tha most wldaly usad and by far 
tha most succassful (Estas at al., 1983). Classification aethods hava 
baan usad In conjunction with teaporal analysis to taka full advantaga of 
saasonal vagatation changas (Rob inova, 1981, Aals at al. , 1981; for nummary 
sea Estas at al., 1983). 

Classification of cropland has rallad on various statistical 

dus taring algor lthas, such as those used in tha Large Area Crop Inventory 
Experiment 

(LACIE) (Lennlngton at al., 1984, Amis at al., 1981). These methods 
generally require tha Identification of training samples or pure pixels, 
l.a. , tha association of Landsat observations with known targets. Each 
training sample consists of a soil and canopy combination with a spectral 
response raprasantatlva of a larger area. This advantaga Is usually not 
present whan classifying natural landscapes, since most targets possess 
unique spectral characteristics. 

The following sections describe several techniques which have evolved 
for analyzing non-unlformly vegetated landscapes. 

Vegetation Indices . Semi-snip Ir leal methods for reducing the amount 
of MSS data have led to the development of the vegetation index (VI). The 
use of the VI for identifying vegetation parameters has received much 
attention In the last few years. Perry et al. (1984) summarize the 
numerous Vi's which have been proposed. A typical VI is a normalized ratio 
of the form, 

VI - * 1R 

“lR + ®VIS 

where a^ and a^g represent the surface albedo in the IR and visible 



‘ ■s ' % ' 



***** tm 


rang**, r«cp«ctlvely . Low VI lndlcacas in vegetation amount, whereas high 
VI indicates elthar high vagatation amount or high productivity (Curran, 
1980). Sallara (1985) discussed th* functional relationship between 
normalized VI and several vegetation parameters, Including th* Leaf Area 
Index (LAI), greenness, and photosynthetic rat*. Tucker at al. (1983) 
correlated a normalized VI (computed from NOAA's AVHRR data) to actual 
biomass obtained from field sampling In a seml-arld region of Senegal, Vest 
Africa. However, much scatter of the data exists. Recently, the effect of 
soil background on Vi's computed from hand-held or airborne radiometer data 
has been evaluated (Elvldge *t al. , 1985, Huet* at al. , 1985). 

While th* functional relationship between th* VI and several 
vegetation parameters has been qualitatively Identified, the quantification 
of these parameters (based on Vi's computed using Landsat data) has 
attained only limited success. Additional research Is warranted In at at 
least two areas. First, the functional relationship between the VI and 

percent canopy cover needs to be clearly established. Secondly, the 

•** 

usefulness of the VI In analyzing natural landscapes needs to be examined. 

Theoretical Models . The problem of resolving a mixed-pixel has been 
approached by assuming that Its total spectral response Is a linear 
combination of the Individual spectral responses of Its components (HcCloy, 
1980; Dozier, 1981; Ungar et al., 1981; Chhlkara, 1984). HcCloy (1980) 
proposed that under conditions of negligible canopy transmission or 
multiple reflection, the response proportions of the various land covers 
will closely approximate the physical proportions of each type of cover. 
HcCloy (1980) suggested that up to four sub-pixel categories be used 


Including three level* of vegetetlon greenness cover end one soil 
beckground cover. 

Unger et el. <1981) reported successful dellneetlon of forest cenopy 
types In Heine using e sleller epproech which they teraed the "Fanning 
algor 1 the". Unger et el (1981) assuaed that the total radiance received 
froa e mixed pixel, I* , can be expressed 

\ - n\ + <1- n ) \ (2) 

where end are the energies received In the X th band froa 
sub-pixel surface eleaents A and B, respectively. The coefficient n Is the 
fractional area occupied by cover type A. The value n was deterained by 
Unger et el (1981) by minimizing the error between observed aud theoretical 
radiances. 

Dozier (1981) also proposed a siailar aethod using two Infrared 
bands for differentiating radiant temperature fields of sub-pixel spatial 
resolution. Corrections for atmospheric effects were Included. 

The theoretical nadir radiance observed by a satellite has been 
derived by several authors (Dave 1980, Otteraan et al. 1979, 1980, and 
Otterman 1973, 1931) for a pure pixel In a background of different 
reflectivity. The nadir radiance Is expressed In terms of three 
components. The direct beam from the object pixel which has been 
attenuated due to atmospheric effects, , Is 

l or X ' r °t •*> <- T )h <3> 

where r ■ reflectivity of object pixel 

G t - global surface lrradlance on object pixel 
t * total vertical optical thickness. 


The total optical thickness Is a sum of optical thicknesses duo to Rayleigh 
scattering by gas molecules, tr, and Mle scattering by aeroools, t 

Assuming a Lambertian surface, the portion of diffuse radiance from 
the surrounding vicinity which Is scattered to the satellite by the 
atmospheric column above the object pixel is 

Hu*" ** G t^ T * r 0 /2 i <" x /cos i )] cos>[t,.P r << ) + t m P m 4 >] 2* sin +# 

where 

a - average reflectivity of the surrounding vicinity 
P R» P H " phase functions associated with Rayleigh and aeroaol 
(Hie) scattering, respectively 
4 ■ zenith reflection angle. 

Finally, the radiance scattered from the direct solar beam back to the 
satellite la written 

L . u Ji-«*pHU+..ce o )}l[T R p R ua° 0 -9o> + <5> 

-fl—, 

where Mo ** cos 6 0 

8 o ■ solar zenith angle 

The total nadir radiance Is the sum of equations (3), (4) and (3). The 
amount of area to be considered in determining the average surrounding 
reflectivity is discussed by Otterman et al. (1979). Otterman et al. (1980) 
reduced the above equations to a simpler fora In the case of an optically 
thin atmosphere. Otterman (1981) also developed an expression for 
reflection from sparsely vegetated soils In terms of a surface cover 
parameter. While accounting for the shadowing effects of vegetation 
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clumps, the formulation assumes that raflactlon from tha top of tha clumps 
is negligible. Thus, thla mathod la not directly applicable to tha 
datarmlnatlon of parcant canopy covar. 

Summary of Lltaratura Review . 

Tha lltaratura review haa lad to tha following conclusions regarding 
tachnlquas available to estimate canopy covar at a subplxal level of 
resolution: 

1) Tachnlquas for evaluating mixed-pixels have focusad primarily on 
distinguishing cropped lands. Most tachnlquas require tha Identification 
of training samples , which prevent thalr application to Inhomogeneous 
canopies. Tha application of these methods to natural landscapes has 
received very limited attentive and needs to be developed further. 

2) While the correlation of the with several vegetation 
parameters , such as blomasci or greenness, has demonstrated definite 
success, no quantitative Information has been found relating the VI to 
percent canopy cover. Further, the use of VI over natural landscapes has 
received little, If any, attention. 

3) The theoretical framework for modeling the observed vegetation 
cover from satellites exists, but has been tested only to a limited 
degree. Ho studies reviewed applied the modeling approach to seml-arld 
systems. Current available models need to be extended to model a mixed 
target pixel. The substantial amount of Information on plant canopy and 
soil reflectances will facilitate testing of a theoretical model. 

Formulation of Percent Canopy Covir Algorithms 

Two approaches for determining the percent canopy cover at a 
sub-pixel scale of resolution have been formulated based on the empirical 


and thaoratlcal relationships reviewed abova. A first approach Involves 
developing an aapirlcal or seal-empirical relationship between ths 
noraallsed vegetation index and the percent canopy cover. A thaoratlcal 
approach Involves adapting the davalopaants of Ottoman (1978, 1981), 
Ottoman at al. (1979, 1980), Dav*» (1980), Ungar at al. (1981), Dozier 
(1981), and McCloy (1980) to a mixed target pixel. Both procedures can be 
verified using ground truth obtained from existing photographs. 

Empirical Approach 

Since the canopy cover is functionally related to other plant and 
soil paraaotars (Rosenberg et al., 1983), It is reasonable to expect good 
correlation between the normalized vegetation Index and percent canopy 
cover. The first approach being considered Is simply to test this 
hypothesis using MSS data and ground truth obtained from existing aerial 
photographs. 

Several variations to the simple correlation scheme are envisioned. 
For Instance, It is possible that the presence of only a few pure or nearly 
pure pixels (either completely bare sofcl or completely green vegetation) 
may be sufficient for providing bounds on the range of the numerical value 
of the VI which occurs within the segment. These bounds would appear at 
the extreme ends of s histogram Indicating the number of occurrences versus 
VI as shown in Figure 1. Once these end points are defined, the relative 
percentage of each cover type can be estimated by methods reviewed earlier. 

The validity of this approach relies on the assumption that there 
exists a clear distinction between soil and canopy reflectances, and that 
the lnhomogeneltles within the two separate classes are smaller than those 
between the soil and canopy. The assumption is more apt to be valid if the 
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Figure 1. Hypothetical Histogram of Number of Pixels 
versus Normalized Vegetation Index 
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segment analyzed la not "large", such that tha variability in diffarant 
■oil and plant types is snail. 

Theoretical Approach 

This procedure treats the direct bean from the target pixel as a 
linear combination of two or three Individual sub-pixel elenents. Equation 
(3) can thus be reformulated for the simple two cover model as 

Lnr X "nl* v + <*-n> Hs < 6 > 

where n ■ the effective portion of a pixel covered with vegetation 
I^v * Radiance emltt^4 from vegetation canopy in band X 
* Radiance emitted from the soil in band X. 

The total nadir radiance will bo the sum of aquations (4), (5), and (6). 

The analyals will be limited to an optically thin atmosphere (clear 
sky) to reduce the Influence of the scattering terms (Equations (4) and 
(5)). Equations (4), (5), and (6) will be evaluated by assuming 
reflectances for the different cover types in each band, as well as the 
clear sky scattering and absorption parameters tr, Pr, t m as 
previously defined. These parameters can be effectively estimated based on 
meteorologlc conditions (Iqbal, 1983). An iterative scheme will be 
developed to minimize the difference between the computed and observed 
radiances, based on successive estimates of n • 

Finally, since there are four MSS bands, it is theoretically possible 
to extend the above procedure to the case of a pixel containing several 
cover types. For example, the presence of trown herbage can significantly 
affect the overall canopy reflectance (Tucker, 1980). To account for this 
third cover type, the expression for the direct beam from a pixel can be 


modified at 


S* ' "KV+ i\i* < 7 > 

where, In addition to tha pravloualy dafinad terms, 

Y - tha affactlva portion of a plxal covarad with brown 
vega tat Ion 

- Radiance emitted froa the brown vegetation In band X* 

Tha use of two bands will be required In this situation. Tha practical 
application of this approach depends on the ability to define unique 
spectral characteristics for each of tha three cover types. This condition 
Is important not only for resolving the data but also to aalntaln linear 
Independence between the two direct beaa equations. 

Selection of Test Sites . 

Significant progress has been achieved toward the selection of 
appropriate test sites for objectives 2 and 3 of the research (see page 1). 
For the testing of the mixed-pixel algorithms, Objective 2, the two major 
criteria are that 1) the region not be fully vegetated at v sub-pixel level 
of resolution and 2) the area must contain Landsat reconnaissance and a 
corresponding aerial photographic survey at approximately the same time. 

The aerial survey must be of sufficient detail that percent vegetation 
cover can be delineated. 

The criteria for selecting a test site for testing Eagleson's 
hypothesis, Objective 3, are more stringent. Test sites must be 
water-limited, naturally vegetated, and at a near climate-climax state. 
Additionally, long term hydrologic and climatologlc records must be 
available. It Is also desirable that climatologlc data during the time of 
the surveys be available. 
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The above constraints hava diraetad tha site-selection process toward 
tha seal-arid regions of southwest United States. Numerous agencies 
located in those regions have been contacted Including the National Park 
Service, the Forest Service, the Soil Conservation Service, the Bureau of 
Land Management, the U.S. Geological Survey, Arizona State Parks, the 
Agricultural Research Service (ARS), and tha University of Arizona. The 
Information provided by the various regional and state offices Is currently 
being examined. However, It appears that the criteria for Objective 2 can 
be easily met. For Objective 3, it Is likely that the eventual site be 
located In the ARS experimental watersheds and/or the National Parks. The 
sites currently under consideration are the Ualnut Gulch Experimental 
Agricultural Watershed, near Tombstone, Arizona, the Big Bend and Guadalupe 
Mountains National Parks, Texas, and the Carlsbad Caverns and Bandeller 
National Parks, New Mexico. These locations are near or in the Sonoran cvc 
Chlhuahuan Deserts. 

The feasibility of obtaining concurrent aerial photographic and 
Landsat data has already been showu to be very encouraging; The selection 
process la being facilitated through the aerial photographic and Landsat 
search services conducted by the EROS Data Center, Sioux Falls, 

South Dakota. Additional sources of aerial photographic surveys which have 
been contacted Include the U.S. Department of Agriculture, Agricultural 
Stabilization and Conservation Service, Salt Lake City, and the U.S. Geo- 
logical Survey, Flagstaff. As of this date, one search request has been 
made through EROS for the Walnut Gulch Hydrologic Experimental Watershed, 
Arizona. Matching of the results for the aerial photographic and Landsat 
searches resulted in one period In October, 1978 In which both data were 


takan within two waaka. Slnca thli parlod is one In which canopy cover la 
not likely to change, both aurveya can be used to teat the aelected canopy 
cover algor itha. Saaple photograph! have been ordered f roa ASCS and EROS 
for tha Walnut Gulch aita. 

Additional aearchea for other areaa are currently underway. The aoat 
effective procedure la firat to aelect the geographic location, then to 
requeat the aerial photographic aearch. If an aerial aurvey of appropriate 
acale can be found for a period later than about 1972, it la very likely 
that correspond lug satellite data will be available In light of the 
frequent coverage of the entire United States by Landaat and NOAA 


satellites. 


FUTURE RESEARCH 


Future research will concentrate on refining the two canopy cover 
formulations outlined above and then applying them to eeveral aelected 
areaf/ In the eouthweatarn United Statea. Once the final altea have been 
aelected, the entire aerial photographic covarage will be purchaaed and 
analyzed for percent canopy cover. The correapondlng Landaat data will 
then be analyzed In accordance with the above propoaed techniques at the 
GSFC. Any refinements to the above procedures will then be made. 

Once a suitable canopy cover procedure is established, It will be 
applied to the appropriate test watershed. Soil hydraulic properties will 
be estimated using Eagleson's hypothesis. Results will be verified using 
actual streamflow data. 

The feasibility of using Landsat data to compute soil hydraulic 
properties will thus be established. A program will then be designed to 
extend this procedure to other natural systems. 
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